The influence of the spin-flip transitions on the 
photon spectrum from ultra-relativistic electrons in 
the field of a crystal. 

Andrei V. Korolfl 0fl, Andrey V. Solov'yovf§ Q, and Walter Greinert Q 

flnstitut fiir Theoretische Physik der Johann Wolfgang Goethe-Universitat, 60054 
Frankfurt am Main, Germany 

^Department of Physics, St. Petersburg State Maritime Technical University, Leninskii 
prospect 101, St. Petersburg 198262, Russia 

§A.F.Ioffe Physical- Technical Institute of the Academy of Sciences of Russia, 
Polytechnicheskaya 26, St. Petersburg 194021, Russia 

Abstract. We investigate the influence of the spin- flip transitions on formation 
of the spectral distribution of the radiation emitted by ultra-relativistic of tens to 
hundreds GeV electrons incident along crystallographic axis in a thin single crystal. 
Both the formalism and the numerical data are presented. The calculated spectra for 
35 to 243 GeV electrons in W crystal are compared with the dependences obtained 
experimentally 



1. Introduction 

The aim of this paper is to carry out quantitative analysis of the role of the spin- 
flip transitions in formation of the total emission spectra by ultra-relativistic electrons 
which traverse single crystal close to its axial direction. This study was stimulated by a 
recent paper |jj where the experimental data on radiative energy loss of e = 35 to 243 
GeV electrons incident on a W single crystal are presented and, to a great extent, are 
interpreted as the first experimental evidence of a significant contribution of spin to the 
radiative spectrum. 

The photon emission by a highly relativistic particle in a static field of crystalline 
axis/plane and the magnetic bremsstrahlung (or synchrotron radiation) are examples 
of the radiative processes occurring in strong fields. It is a well-established fact (e.g. 
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|2|, that quantum effects start manifesting themselves in the photon spectrum when a 
dimensionless quantum parameter \ exceeds the value of ~ 0.1. The quantum parameter 
can be written as x = 1 H/Hq (for magnetic bremsstrahlung) or x — lE/Eo (for the 
radiation process in the external electric field), where 7 = e/mc 2 is the relativistic factor 
of projectile, H/E is the strength of the magnetic/electric field and Hq = Eq = m 2 c 3 /eh 
is the critical field. 

In this paper we focus attention on the influence of the spin-related effects on 
the formation of the spectral distribution of unpolarized photons emitted by tens to 
hundreds GeV electrons incident under small angle $0 to crystallographic axes of a thin 
single crystal. 

Theoretical approaches suitable for the analysis of the influence of the spin-related 
effects on the radiative spectrum of ultra-relativistic particles in a strong external 
field include (a) the quasi-classical method based on the operator representation of 
the projectile's wave function ||, || (see also |J), (b) the approach based on the 
correspondence principle and utilizing the equivalent photon method (see also ||), 
(c) the full quantum-electrodynamical treatment of a kinetic equation for an electron 
moving in the field of a crystal axis . 

In section |2| we present general formulae which allow to estimate quantitatively the 
contribution of the spin-flip transitions to the total radiative spectrum. This part of the 
paper, which is included to make the reading easier, contains all the essential expressions 
which were used to carry out the numerical analysis. They represent the quasi-classical 
formalism described in great detail in ||. 

To say it very clearly: (a) All the formulae are from Q except one (equation (|9|) 
below) which could not be found there. This refers to the explicit expression for 
the differential intensity with spin-flip, (b) In application to the radiation emission 
in crystals the formulae contain the standard CFA (constant field approximation) and 
correction terms (for details see section H|). (c) When carrying out the numerical analysis 
we used both methods, 'pure CFA' and 'CFA + corrections'. 

The numerical results are presented, discussed and compared with the experimental 
data Jj]] in section || A summary of the results and possible directions for further 
improvement of the theoretical model are given in section [|. 



2. The formalism 



Within the framework of the quasi-classical approach the spectral-angular distribution 
of the energy radiated by an ultra-relativistic particle (of a spin s = 1/2) moving in an 
external static field is given by the following formulae 
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The notations used are as follows: a is the fine structure constant, u, k and e are, 
respectively, the photon energy, wave vector and polarization vector, n = ck/w denotes 
the direction of the emission, c is the light velocity, df2 n is the solid angle associated 
with n. The initial and final energies of the projectile are denoted by e and e' = e — hu, 
the quantities £ and £' stand for the polarization vectors of the initial and final states 
respectively, a is the vector constructed from the Pauli matrices <jj (j = x,y,z). The 
quantity to' is defined as uj' = uj {e/e') 

The remarkable feature of the quasi-classical formulae (HHJ) is that they combine 
the classical description of the projectile's dynamics and the specific quantum effects. 
Indeed, the integrand in ([!]) is time- dependent through the functions p(t) and r(t) 
which represent the classically defined momentum and the position vector of the 
projectile moving in an external field. Therefore the emitted radiation is described 
in terms of the trajectory defined by p(i) and r(i), as in classical electrodynamics; it is 
uniquely parametrized by the initial energy e and and the initial values po = p(0) and 
r = r(0). On the other hand, these expressions account for two important types 
of quantum corrections: (a) the radiative recoil (i.e. the change in the projectile 
energy /moment a due to photon emission) which formally enters the integrand in ([!]) 
through the 'corrected' photon frequency uj', and (b) the dependence of the intensity of 
radiation on the spin quantum numbers of the projectile. The latter enter the formulae 
through the two-component polarization density matrices (1 + £ • a)/2 and (1 + £' • <?)/2 
(see©). 

The quantum corrections related to the projectile's dynamics in the external field 
are neglected. This is fully justified provided the particle can be considered as an 
ultra-relativistic in both the initial and the final states. Hence, it is implied that the 
relativistic factors satisfy the strong inequalities 7 = e/mc 2 3> 1, 7' = e'/mc 2 3> 1. If 
these conditions are met then, as it is proved in 0, Q (see also the recent paper 0), the 
expressions tf^fy can be used to evaluate the characteristics of radiation emitted by a 
projectile moving in an arbitrary static field. 

Once the functions p(£) and r(t) are known, then the these formulae allow one 
to evaluate the most detailed characteristics of the emitted radiation. These include 
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the dependence of dE not only on the photon energy and the emission angle but the 
polarizational properties as well. In particular, one can analyze the role of the spin 
degree of freedom in the radiative process. 

To obtain expressions describing the radiative spectra of a collimated bunch of 
unpolarized electrons passing through a crystal close to the axial direction one has to 
carry out the following transformation of the right-hand side of ([I]): (a) summation 
over the photon polarizations, (b) averaging over the electron polarizations (the cases 
£' = ±£ have to be considered separately), (c) integration over the emission angles, and 
(d) averaging over the fast oscillatory transverse motion of the electrons due to the action 
of the axial field. The transformations (a)-(c) are carried out using standard methods 
(see, e.g., H, f|). An adequate approach for averaging over the transverse motion is based 
on the statistical description of the particle distribution in the phase space (|5|) |10fl). In a 
thin crystal, where one can neglect the influence of multiple scattering from the crystal 
nuclei and electrons, the differential probability, dw, for a projectile with (transverse) 
coordinate p and velocity can be written as follows 

dw(p,v ± ) = ^F(p,$ ), (6) 
s 

where s is the area corresponding to one axis, and $o is the incident angle with respect 
to the axis (see the illustrative figure 1).' The distribution function F (p, $ ) is defined 




Figure 1. Schematic representation of the passage of an ultra-relativistic particle 
through an oriented crystal. Under the action of the crystalline field U (the filled 
circles denote the positions of atomic nuclei arranged into the crystallographic axes, 
- atomic strings) the trajectory r(t) of the particle varies along the direction of the 
vector v , which is the mean velocity. The picture corresponds to the case when 
incident angle #q exceeds the Lindhard's angle d^. The dashed ellipse denotes the 
area, s — irp^, allocated for one axis. The field of each atomic string is supposed to be 
axially symmetric within s: U = U{p) where p is the distance from the axis. 

by the initial conditions at the entrance. In the case of an axially symmetric potential 



5 



U(p) of the atomic string the function F (p, $ ) = F (p, i? ) is given by: 

F(p,# )= J ^^9(e ±0 -U(p)), (7) 

where s (ej_o) — / d 2 p 6 (e±o — U(p)) is the allowed area for the transverse motion in a 
unit cell, 6(x) is the Heaviside function: 9(x) = 1 when x > 0, and 9(x) = if otherwise. 
The quantity e±o stands for the transverse energy e±o = £$l/2 + U(po) which, for given 
t?o and the entrance coordinate po, is the constant of motion in a thin crystal. Note that 
for the incident angles larger than the critical Lindhard-angle, i?£ = (2Uq/s) 1 / 2 (where 
Uq is the depth of the potential well), when e±o > Uq, the right-hand side of ([?]) reduces 
to unity. 

Carrying out all the operations mentioned above one obtains the spectral 
distributions of the emitted energy (per path dl = c dt) in the following form: 



where the brackets designate the averaging over the phase space, and the bars indicate 
that the averaging over the initial polarizations of the projectile is carried out. The 
integration is carried out over the area allocated for one axis, s = irp 2 , where p s is half 
the distance between two neighbouring axes of the same type. 

The superscript '(t)' specifies the type of the spectral distribution: (t) = (±) 
correspond to the spectra with, (-), and without, (+), the spin-flip transitions, and 
(t) = (tot) designates the total spectrum which is the sum of the (+)- and (-)-terms. 

The spectral distribution appearing in the integrand of (|D reads 
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where A c is the Compton wavelength of the electron. For the cases '(t)=(tot)' and 
'(t) = (-)' the coefficients and (5^ are given in terms of the fractional photon energy 
fi — tku/e, < p, < 1: 

^ (tot) = -/i, C^ = -^— v (3^ = 1-11 + -J- , /3H = -2. (10) 

o(l — p.) 1 — p. 

For small incident angles one can use the continuous model for the potential of the 
axis. Then, assuming this potential to be axially symmetric, solving (approximately) 
the equation of motion v(t) = (dU(p)/dp) using on the right-hand side the relation 
r(i) = ro + vot for the position vector (vo is the the mean velocity of the particle), and, 
finally, using the obtained expressions in eqs. (l)-(5), one finds that the functions 3>(t) 
and A(t) can be represented as follows 
52 „ ^ „l 
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where d v = Uo/mc 2 , W(q) is the atomic formfactor, ut is the amplitude of thermal 
vibrations of the atoms, and go = q • vo. The integration is carried out over the vectors 
q and q' lying in the plane perpendicular to the axis. The details of derivation of eq. 
(|) in the case of the total spectrum can be found in || . 

Let us note here that the radiative transitions accompanied by the spin-flip process 
manifest themselves for large energies of the emitted photon, hu ~ e. Indeed, it follows 
from (|9|-|lOD that dE(~~) oc /i 3 /(l — /i). The latter quantity is negligibly small for /i < 1 
and increases as \i ~ 1. It is known (see, e.g. 0, that the emission of such energetic 
quanta becomes noticeable in the radiative spectrum when the quantum parameter Xs 
becomes large, Xs > !• Hence, the magnitude of Xs defines the effective range of the 
emitted photon energies. In the case of an ultra-relativistic particle passing through a 
crystal close to the direction of a crystallographic axis the parameter Xs is conveniently 
represented in the following form: 

A c Uo , s 

Xs = l 2- (13) 

a s mc z 

The quantity a s represents the effective screening radius of the potential created by a 
string of the crystal atoms. For the potential of (111) axis in a tungsten crystal, where 
U ~ 400 eV, a s ~ 0.2 A(see, e.g. [§J), it follows that Xs = 1 is achieved for the incident 
electron energy 35 GeV. Hence, one can expect that the effects due to the spin-flip 
transitions will manifests themselves in the total radiative spectrum of an e > 35 GeV 
electron moving at small angle $o with respect to (111) direction in a W crystal. 

^From the computational viewpoint the direct calculation of the functions ( 1 1 -|T^) 



and, correspondingly, the spectra (|9|) is a non-trivial problem even if one utilizes simple 
analytic forms for the axial potential such as the Moliere approximation or the Lindhard 
potential [1C]. Considerable simplification of the formulae is achieved within a so- 
called constant field approximation (CFA) |J. The latter implies that the change in 
the particle's radial distance, Ap = v±r c = c&qt c , occurring during the time interval 
T c — lc/c (here l c ~ cy 2 /u' is a coherence length), is smaller than the typical scale a s 
within which the axial field changes noticeably. Hence, provided the condition l c i9o < a s 
is met, one can expand the integrands in ( pjjT2] ) in powers of g r ~ c&Qr/a s < 1. If, in 
addition, the incident angle $o is small compared with then the terms proportional 
to (f?o/^i)) 2 can be treated as small corrections. Within this approximation the functions 
( P^[fT2] ) are expressed in terms of local characteristics of the trajectory. The final result, 
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valid for both the total spectrum and for its part due to the spin-flip, reads 

where Ai(.) and Ai'(.) are the Airy function and its derivative respectively (see, e.g., 



111). Other variables and functions are defined as follows: 
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By expressing Ai(z) and Ai'(,z) in terms of the modified Bessel functions Ki/ 3 (x) 
and K 2 / 3 {x) with x = 2z 3 ^ 2 /3 (see, e.g., |nj) we reproduce the formula presented in §17 
of[|. 

The first two terms on the right-hand side of (fH|) are independent of the incident 
angle 9q and represent the quasi-classical spectral distribution of the radiated energy in 
the magnetic bremsstrahlung limit (see, e.g. [@, 0). 

The last three terms on the right-hand side of ([14]) represent the correction to the 
CFA due to the change in the field strength on the scale of the coherence length, l c . 
These terms appear as a result of the transformation of a more general expression (|9|). 
The range of validity of ([TJ]) is subject to two conditions: (a) $ and (b) Z c i? < <V 



3. Numerical results 



In this section we present the results of calculations of the spectral intensity of radiation 
emitted by ultra-relativistic electrons in the energy range e = 35—243 GeV incident with 
$o = 0.3 mrad from (111) axis in a W crystal of a thickness 0.2 mm. The parameters 
correspond to those used in a recent experiment |IJ where the main emphasis was put 
on the investigation of the role of the spin-flip process in forming the total radiative 
spectrum. 
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The calculation of the spectral intensity was based on the formulae ( p^ - |20| ). The 
averaging over the projectile trajectories was carried out according to the rule (|8]) where 
we used F (p, $o) — 1 which is justified by the relation i9 Q > Indeed, for the potential 
of a (111) axis the depth of the potential well is ~ 1000 eV and, therefore, the Lindhard 
critical angle d L = 0.09 ... 0.2 mrad for e = 35 - 243 GeV. 

We used a continuous potential approximation for modelling the field of an atomic 
string. The potential of an individual atom was calculated by using the radial 
wavefunctions within the non-relativistic Hartree-Fock approximation. For the distances 
p = [0,p s ] (where p s = L c /v6 is the half-distance between two (111) axes, L c = 3.165 
A is the lattice constant of a W crystal ||12|| ) the potential of the individual string was 
corrected by accounting for the potentials created by the nearest six axes. Only axially- 
symmetric part of the resulting field was used. The dependence of the field strength 
dU(p)/dp on the distance p from the axis are plotted in figure |[ The screening distances, 
a s , defined as the root of the equation dU(a s )/dp = e -1 (dU/dp) max , are equal to 0.125 
A for the crystal temperature T = 70 K and to 0.191 A for T = 293 K. 

As functions of p the spectral intensities dE^/hdu dl depend strongly on the 
magnitude of the local quantum parameter x(p) defined by eq. (0). Namely, the 
value of x{p) defines the effective range of the emitted photon energies. The formal 
analysis of this statement is as follows. The right-hand side of ( |I~4"D is expressed in 
terms of the Airy function and its derivative. These functions satisfy the conditions 
Ai(z), |Ai'(z)| = 0.1 ... 1 for z < 1, and decrease rapidly (exponentially) as z goes beyond 
1. The same is valid for the integral term, / z °°Ai(£) d£ (see, e.g., JIT]]). Hence, in the p- 
range where x{p) ^ 1; the spectral distributions increases sharply within the interval 
of low energies of the emitted photon, Tiuj = . . . huj max with Tiu max « x(p) £ ^ £ > 
and decrease rapidly as Tiuj/e > x(p)- If x(p) ^ 1 then the condition z ~ 1 produces 
(see (|i~5l) ) hu max /e ~ 1 which makes possible the emission of energetic quanta. In this 
case the total spectral distribution d£'( tot )/^da;d/ is (roughly) flat within the interval 
huj/e = 0...1 0. The contribution of the spin-flip term becomes noticeable (and 
comparable in magnitude with the total spectrum) at the hard end of the spectrum 
whereas for hu/e 1 it is negligibly small. The latter result immediately follows if 
one compares the coefficients C^ *-* and (see (|10|)): for p 1 the following strong 
inequality is valid C^/C^ w p 2 < 1. 

Figures |3](a-b) represent the dependences x{p) f° r t wo temperatures of the crystal 
and for different energies of the projectile (as specified). For the sake of comparison 
the values of the quantum parameter Xs (see fljjD) are plotted as well. The maximum 
distance p s is defined by the crystal symmetry and the chosen crystallographic direction 
and is independent on the crystal temperature. For a (111) axis in a W crystal p s = 1.292 
A. As mentioned above the shape of the spectral distribution ( |I~4"D is quite sensitive to 
the magnitude of the local quantum parameter x(p)- The range pi < p < p2 (where 
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Figure 2. The field strength, dU(p)/dp, of a (111) axis in W as a function of the 
distance from the axis for two different temperatures as indicated. 



p 12 are the roots of the equation x(p) — 1) establishes, for given temperature and 
energy, the distances from the axis which contribute effectively to the emission of high 
energy photons {hu ~ e) and, consequently, define the fraction of spectrum due to 
the spin-flip transitions. For fixed temperature the values of pi 2 do not depend on e, 
while the interval [pi,p2] increases with T. Thus, the intervals of the radial distances 
which contribute to the spin-flip part of the total spectrum ranges from rs [0.05,0.18] 
A for T = 70 K to « [0.0, 0.44] A for T = 293 K. Comparing these values with the 
range p < p s = 1.292, one can anticipate that the radiative spectrum averaged over the 
whole interval of the radial distances (see (|j)) is influenced, to a great extent, by the 
contribution of the regions where x{p) < 1- As a result, even for high energies of the 
projectile electron the relative yield of soft photons is larger than that of %u ~ e. 

Keeping in mind the large variation in the values of the local quantum parameter 
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Figure 3. The local quantum parameter x(p) ( see (PlO versus the distance from the 
axis p for the field of (f If ) axis in W calculated for T = 70 (figure (a)) and for T = 293 
K (figure (b)), and for different incident energies: e — 35 GeV (solid lines), e = 70 GeV 
(dotted lines), e — lfO GeV (dashed lines), e = f49 GeV (long-dashed lines), e = 243 
GeV (dot-dashed lines). The horizontal lines mark the corresponding magnitudes of 

Xs- 



within the interval p < p s of the accessible distances (x(p) is of the order 10 -3 for 
p ~ p s , - this range is not presented in figures §(a-b)), it becomes clear that to obtain 
realistic estimates of the spectral distribution of radiation one has to carry out the 
proper averaging procedure. 

The results of our calculations of the radiative spectra are presented in figures |] 
(a-e). The energies of electrons as well as the incident geometry (as specified in the 
caption) correspond to the experimental conditions |IJ. In addition to the information 
presented in the caption let us note the following. The averaging procedure was carried 
out according to (|8|) where we used F (p,-d ) = 1 which is adequate for the over-barrier 
motion in a thin crystal. Thus, we have disregarded the evolution of the distribution 
function due to the effects of multiple scattering and the radiation damping. All curves 
in figures [| (a-e) were obtained within the CFA, i.e. by using the first two terms in 
the curly brackets on the right-hand side of (0). The dependences represented by open 
circles were calculated (at T = 293 K) using a more elaborated method. For the photon 
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frequencies satisfying the condition l c $ < a s we used flT4] ) including the correction terms 
proportional to To calculate the total radiative spectrum in the region of lower 

photon energies we made use of the asymptotic expressions (17.11-12) from 0] which 
are valid, within logarithmic accuracy, for Xs 3> 1 and p > Xs (a s /p s ) 3 ^ 2 - 

For the sake of comparison in figure |](f) we plotted the spectral distributions 
calculated for a e = 243 GeV electron and for the incident geometry as described 
above but without averaging over the radial distances p. The curves corresponding to 
T = 70 and 293 K were obtained by accounting for the first two terms in QUip where the 
constant values of the quantum parameter x = Xs (see figure H) were used. It is exactly 
this approximation which was used, as one understands, in JE] when interpreting the 
experimental data. 

First we discuss the differences between the averaged and non-averaged intensities. 
There are three features which clearly distinguish the curves in figure |] (f) from their 
analogies in figure § (e). These are: (a) the absolute magnitude of the intensities, (b) 
the shapes of the spectral distributions, and (c) the differences between the T = 70 K 
and T = 243 K spectra. 

Although striking, these differences can be easily understood. To illustrate this let 
us analyze the behaviour of the intensities due to the spin-flip transitions. The analysis 
of the discrepancies in the total spectral distributions can be carried out using similar 
arguments. 

It is easily checked that for large values of Xs the corresponding non-averaged 
spectral intensity due to the spin-flip transition attains its maximum at the photon 
energies satisfying the condition u — p/(l — p) ~ Xs- A quantitative estimate of the 
intensity in the vicinity of the maximum reads 

X-Ai'(l), (21) 



h da; dl I x=Xs 3A C 7 



where Ai'(l) « -0.16. 

To estimate the magnitude of the averaged intensity in the photon energy range 
u m Xs one first notices that the principal contribution comes from the spatial region 
where the local quantum parameter x{p) is greater than Xs (see figure ^). Denoting the 
roots of the equation \{p) — Xs through pi 2 , and taking into account that p\ -C p\ one 
obtains the following estimate: 



/ d#H \ ^ a fP2 2pdpAi'(z) n 



\ hduj dl J u ^ Xs 3A C 7 J P1 pj z 3A C 7 

where (x) is the average value of the quantum parameter over the region pi < p < p 2 . 
For the temperatures 70 and 293 K (x) equals to 21.1 and 12.3, respectively. The 
quantity (z), defined as (z) = (x s / (x)) 2//3 > equals 0.65 for T = 70 K and 0.71 for T = 293 
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(a) £=35 GeV 
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Figure 4. In figures (a-e): the averaged spectral distributions of radiation (see (||)) 
calculated for different energies of electrons (as indicated) aligned 0.3 mrad from 
the (111) axis in W. Thick lines refer to the crystal temperature 293 K, thin lines 
correspond to T — 70 K. Solid lines and open circles represent the dependences of 
the total intensity, dotted lines stand for the contribution of the spin- flip transitions, 
dashed lines reflect (dE( + ) /%dw dl). See also explanations in the text. Figure (f) 
shows the spectral distributions of radiation from a e — 243 GeV electron in W at 
T = 70 and 243 K calculated without the averaging over the p-interval and for fixed 
values of the quantum parameter, \ — Xs- The legend is as in figures (a-e). 
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K. One easily checks (see, e.g., [[11]]) that for both of these values (z)^ 2 Ai' ((z)) « —0.167. 
The parameter r\ stands for the ratio p\j ' p 2 s yielding the values rj « 0.012 for T = 70 K 
and rj fa 0.019 for T = 293 K. The factor Xeff = V (x) represents by itself the effective 
quantum parameter which defines the absolute magnitude of the averaged intensity in 
the range of the photon frequencies hu ~ e. 

The right-hand sides of eqs. ([H]) and (p2|) , reproducing with a reasonable accuracy 
the maximum values of the spectral intensities due to the spin-flip transition (the dotted 
curves in figures | (e) and (f)), allows one to explain the discrepancies between the 
averaged and non-averaged spectra. The excess of the non-averaged intensities over 
the averaged ones is reproduced by the ratio Xs/Xeff which equals 44 at T = 70 K 
and 31 at T = 293 K which are close to the exact values which follow from figures 
f| (e) and (f): 57 and 43, respectively. The different behaviour of the two types of 
the spectral distributions also follows from the above written formulae. Indeed, the 
ratio Xs (T = 70) j\s (T = 293) = 1.7 explains the difference between the thick and 
thin dotted lines in figures £| (e), whereas the ratio Xeff (T = 70) /Xeff (T = 293) ~ 1.1 
illustrates that the averaged intensities are much less sensitive to the change in the 
crystal temperature. 

The behaviour of (dEK tot )/7ida; dl) , represented by solid lines in figures (| (a-e), 
illustrates the known tendencies of the total spectral distributions obtained within the 
CFA ||. Namely, for comparatively low incident energies, when the local quantum 
parameter is small in a wide interval of the distances (see figure 0) the spectrum has 
maximum at hu <C e and decreases rapidly with the photon energy. As e increases, 
the range where x{p) ^ 1 increases as well giving rise to the yield of the photons with 
hu ~ e. The magnitude of the intensity in the low-o; region decreases and the spectrum 
becomes flatter over the whole range of photon energies. 

Similar features characterize the total spectral distributions obtained from (j|) and 
(H]) including the correction terms proportional to $0/1?^ ~ 0.16, - the open circles in 
figures |3] (a-e). The contribution of the correction terms is negligibly small for hu ~ e 
when the coherence length is small enough so that the characteristics of the spectrum are 
defined by the instant value of the particle's acceleration. Their influence becomes visible 
in the low-energy part of the spectrum. As it is seen, the inclusion of the correction 
terms decreases the magnitude of the spectrum for hu ^ e which leads to a further 
flattening of the distribution. 

Commenting on the spectral intensities due to the spin-flip transitions we note the 
following. For all incident energies the spectra become pronounced close to the tip- 
region. This feature, which is known for both the synchrotron radiation [0, |3], |J and for 
the radiation by electrons in the field of crystals |], |, 0, 0, can be easily traced from 
the general formulae (|9|-|10D. Indeed, the spin-flip intensity is proportional to the factor 
/x 3 /(l — fi) which is small for /i < 1 and increases as fi ~ 1. For fi — > 1 the spectrum 
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dE(~}/h du> dl exponentially decreases due to the highly oscillatory exponential factor 
in the integrand in @. The contribution to the averaged spectrum (§) comes from the 
spatial region where x{p) ^ 1- The position of the maximum of (dE^ /hdu dl) can be 
estimated as /i ~ 1 — (x) -1 , see (22). The latter formula, as it is mentioned above, allows 
to estimate the magnitude of the maximum intensity, which is almost independent of 
7 for sufficiently high energy of the electron. Indeed, apart from the explicit factor in 
the denominator, the dependence on 7 enters the right-hand side of (|22| ) through (x) 
and (z). Since the ratio (x)/l is independent of the relativistic factor (see (0) and 
(|16l)), the dependence in question is concentrated in (2;) 1//2 Ai' ((z)) which is a slowly 
increasing function of 7. This feature is clearly seen in the figures: the maximum of 
(dE(~^ /h du) dl) increases by a factor of 2 with the energy variation from 35 to 243 GeV. 

Let us mention that the spin-flip spectral distributions were obtained within the 
CFA, omitting the correction terms on the right-hand side of (|T4"D. These terms 
contribute noticeably to the region /i C 1 but are negligibly small closer to the hard-end 
of the spectrum. 

Finally, let us present the numerical data on the effective radiation length, L e ff. 
This quantity, which is defined as 

-, (23) 




defines the interval within which a particle loses its energy due to the radiation. 

The values of L~j^ calculated using the averaged CFA spectra with the correction 
terms included (open circles in the figures), and the ratios r = L'^/L]^, where 
L&h ~ ^Z 2 otr\n In (l83Z~ 1 ^ 3 ^j (n is the volume density of the crystal atoms, r$ is 
the electron classical radius, Z is the atomic number) is the radiation length due to 
bremsstrahlung in an amorphous medium, are summarized in table |l[ The values of r 
are in good agreement with the measured data |1[ except for the case e = 35 GeV where 
our result is ~ 10% higher. 



Table 1. The values of L}, (see (p3|)) and r = L e h/L B ^j calculated for different 
energies of electron. 

s (GeV) 35 70 110 149 243 
L~} } (cm- 1 ) 21.50 23.89 25.64 25.77 25.86 
r 7.10 7.88 8.46 8.51 8.54 



In figure [5] we present the results of our calculation together with the experimentally 
measured spectra [|T|. The data refer to the scaled spectral intensity, i.e. (dE^/hduj dl) 
divided by L~j^. The crystal temperature is chosen to be 293 K. It is this T-value which 
had been used in the experiment |]T| as one can deduce from the paper basing on the 
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Figure 5. Scaled intensities for different incident energies: (a) e = 35 GeV, (b) 
e = 70 GeV, (c) e = 110 GeV, (d) e = 149 GeV, (e) e = 243 GeV. Thick solid 
lines correspond to the total intensities within CFA. Thick dot-dashed lines represent 
the total intensities and thick dashed lines stand for the intensities without spin-flip 
both calculated within the approximation 'CFA + correction terms'. Thin solid lines 
denote the distributions calculated without the averaging over the p-interval and for 
fixed values of the quantum parameter, \ = \s- All curves correspond to T = 293 K. 
Full circles are the experimental data (lj. 
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cited values of the parameter Xs- Let us note that we failed to reproduce the curves 
presented in the cited paper and which, as it follows from the text, correspond to the 
CFA approximation with constant values x = Xs without averaging over the distances 
p (full lines in figure 1 in ]![). In figure [| thin lines correspond to the calculation of 
dE( tot ) /hdu dl according to (13) but omitting the terms proportional to ^q/'&I and using 
the values Xs (as indicated in figure 0(b)) instead of the local parameter x(p). Although 
it seems like the same approximation is used the shapes of these curves noticeably differs 
from those presented in the cited paper. In any case this approximation, in our mind, 
is not adequate to describe the radiative spectra in the crystal (see the discussion in 
connection with figure |J), so that the thin solid curves are plotted for the illustrative 
purposes only. 

'The scattering of the experimental points indicate the uncertainties' Taking 
this into account we may state that for jj, > 0.7 both thick curves, which describe the 
total intensities, match the experimental data. Note that it is exactly the region where 
the spin-flip transitions contribute to the total spectrum (see figure |4j). In the mid- 
range of the //-values the computed intensities underestimate the experimental ones, 
and for /i 1 theoretical results noticeably overestimate the measured spectra. All 
these features are typical, as it is seen from the figures, for the whole range of the 
incident energies, and, in general, the calculated curves are skewed to the left wing of 
the spectrum as compared to the measured data. 

For the sake of comparison we plotted also the scaled intensities of the radiation 
excluding the contribution of the spin-flip transitions (dashed lines). These were 



obtained as L e fj y{dE^ ot y h duj dl) — (dE^/hdu dl)j with both terms in the brackets 
calculated with the help of (§) and (H), and retaining the correction terms in the latter 
formula. Taking into account the experimental uncertainties one can observe that these 
curves taken alone agree fairly well with the measured data on the right wing of the 
spectrum. 

Thus, on the basis of these calculations we cannot confirm in full the statement 
made in |l[ that the region jj, ~ 1 of the radiative spectrum by ultra-relativistic electrons 
moving under a small incident angle with respect to the crystallographic axis is formed 
mainly due to the spin-flip transition of the projectile. 



4. Concluding remarks 

We have calculated the spectral distributions of radiation formed during the passage of 
multi-GeV unpolarized electrons along (111) axis in W crystal. The consideration was 
carried out within the framework of the quasi-classical approach || which adequately 
describes the radiative process of ultra-relativistic particles in strong external fields. 
Special attention was paid to the investigation of the contribution of the spin-flip 
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transitions of the projectile to the total spectrum. The results of our study can be 
summarized as follows. The role of the spin-flip transitions in forming the hard-photon 
end, Too ~ e, of the radiative spectrum is quite noticeable. Definitely, it is much 
larger than that predicted in |9j where a full quantum-electrodynamical approach was 
developed. However, we also can hardly agree with the statement made in |IJ that it is 
the spin-flip transitions which do modify dramatically the shape of spectral distributions. 
The numerical results presented in this paper agree well with the experimental data in 
the range of photon energies %uo ~ (0.7. . . 1.0)e, i.e. exactly where the contribution of 
the spin-flip transitions becomes pronounced. This fact demonstrates the applicability 
of the standard CFA model (combined with a proper procedure for averaging over the 
phase volume of a projectile) for the description of the hard-end of the spectrum for over- 
barrier particles. For lower photon energies, when the coherence length becomes larger 
than the typical scale of variation of the potential of a crystallographic axis/plane, it is 
necessary to go beyond the CFA scheme. In the present paper we utilized the approach 
suggested in || which allows to incorporate the corrections proportional to ($ /$ v ) 2 . 
Our data shows that even when this parameter is small (in the case considered in this 
paper (fto/fly) 2 = 0.16) the correction terms reduce noticeably the magnitude of the 
spectral intensity in the range of small photon energies. This approach becomes less 
adequate as liw/e < 1. In this case one has to use the formulae (0-0) which include 
the deviations from the CFA in the most general way. 

The theoretical results can be improved further by taking into account three 
important effects which are intrinsic for the radiative scattering of an ultra-relativistic 
particle in the field of a crystal. The first one is the analogue of the Landau- 
Pomeranchuk-Migdal effect ]T3 |. The latter consists in a strong suppression of the yield 
of the soft-photon bremsstrahlung emitted by an ultra-relativistic particle in medium. 
This is a result of a number of correlated small-angle collisions of the projectile with 
the lattice atoms on the scale of the coherence length. Similar effect occurs when 
the particle moves close to an atomic string in a crystal and under the small angle 



$ from the axis fl4}| . The account for this mechanism will reduce the magnitude of 
the spectral intensity for low photon energies. Another phenomenon to be accounted 
for is the cascade process of the photon emission which is effectively generated when 
a multi-GeV electron/positron enters the crystal at small ($0 < angles to the 



crystallographic axis/plane [||, I5|. The third effect we would like to mention is the 
influence of the radiative loss on the shape of the emission spectrum. In the present 
calculations we disregarded the change in the particle's energy on the scale of the crystal 
thickness L. However, as soon as L becomes comparable with the effective radiation 
length the formalism must include the solution of the proper kinetic equation (|5|, £|). 
This consideration is fully appropriate in the case considered in the present paper with 
regard to the experiment [|TJ. Indeed, the experimental setup included the 0.2 mm 



18 



thick W crystal. Comparing this length with the values of L e ff (see table [I]) one finds 
L/L e ff ~ 0.5, so that the electron looses half of its initial energy. Although it is 
clear, qualitatively, how the radiative loss influences the radiative spectrum (the latter 
is shifted towards the low-energy part), the rigorous quantitative treatment is to be 
performed for each specific set of the parameters: the incident energy and the geometry, 
and the type of the crystal. 

Finally, let us mention the phenomenon which was left beyond the scope of this 
paper and which, as far as we know, has not been studied neither theoretically nor 
experimentally. It is well-known (e.g. |§) that the magnetic bremsstrahlung leads to 
the radiative self-polarization of the electron. Formally, this is due to the fact that the 
differential probability contains the term proportional to £ • H (here H is the vector of 
the external magnetic field) which depends on the angle between £ and H. A similar 
effect was studied in [TJJ in connection with the axial channeling of electrons/positrons. 
If one takes into account the radiative loss (which is the case when x > 1 an d 7 S> 1) 
then there appears a strong interaction of the spin with the field of radiation, H r 
0j . This results in (a) a noticeable change of the typical time r of the radiative self- 
polarization, and (b) in radiative damping which leads to the decrease of |£|. Later these 
phenomena were considered in connection with the planar channeling of ultra-relativistic 
electrons/positrons in slightly bent crystals fl?| , [378] , |19|| . It was demonstrated that there 
appear additional terms, due to the centrifugal force in the bent channel, which influence 
the spin-precession, and, as a result, the time of the self-polarization. We note here that 
the planar channeling is more preferable than the axial one if studying the effect of 
radiative self-polarization. The reason is that the terms proportional to £H (H is the 
effective magnetic field acting on the particle in the channel) are less influenced by the 
averaging over the phase distribution in the planar case. 

In our mind it will be very interesting to carry out the analysis (theoretically and, 
if feasible, experimentally as well) of the self-polarization of a bunch of ultra-relativistic 
particles undergoing channeling in a crystal whose channels are periodically bent. In 
this case there arises very intensive radiation, additional to the well-known channeling 
radiation, due to the periodicity of the trajectory of the projectile p0 . 
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